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' We report on numerical studies of the NLO corrections to exclusive meson electropro- 

, duction, both in collider and fixed-target kinematics. Corrections are found to be huge 

' at small xb and sizeable at intermediate or large xb- 

^ ■ 1 Motivation and general framework 



Generalized parton distributions (GPDs) are a versatile tool to quantify important aspects 
of hadron structure in QCD. They contain unique information on the spatial distribution of 
partons [1] and on the orbital angular momentum they carry in the proton [5]. The theo- 
retically cleanest process where GPDs can be studied is deeply virtual Compton scattering 
(similar to inclusive DIS, which plays a dominant role in constraining the usual parton den- 
Oh| sities) . Hard exclusive meson production is harder to describe quantitatively, but it provides 

opportunities to obtain important complementary constraints. In particular, vector meson 
production is more directly sensitive to the gluon distributions, which enter the Compton 
amplitude only at next-to-leading (NLO) order in a^. Together with a wealth of high-quality 
I data [3] , this warrants efforts to bring meson production under theoretical control as much 

as possible. 

I In the present contribution [4] we investigate exclusive p production {'j*p pp) using 

collinear factorization, which is applicable in the limit of large photon virtuality at fixed 
CO ' Bjorken variable xb and fixed invariant momentum transfer t to the proton [S]. In practical 

i terms, this means that the description is restricted to sufficiently large but can be used 

' for both small and large x b , thus providing a common framework for analyzing both collider 

. and fixed-target data. The process amplitude can then be expressed in terms of GPDs for 

the proton, the qq distribution amplitude for the p, and hard-scattering kernels. The kernels 
are known to NLO, i.e. to order 

The requirement of "sufficiently large" is demanding for meson production. Contri- 
butions that are formally suppressed by powers of cannot be calculated in a completely 
' systematic way, but the estimates [II El [9] agree that for of several GeV^ the effect of the 

transverse quark momentum inside the meson cannot be neglected in the hard-scattering 
subprocess, as it is done in the collinear approximation. This effect can be incorporated 
in the modified hard-scattering picture [71 [H|j in color dipole models [5], or in the MRT 
approach |10| . Unfortunately, the calculation of corrections remains not only a technical 
but even a conceptual challenge in these approaches, so that the perturbative stability of 
their results cannot be investigated at present. One strategy in this situation is to study the 
NLO corrections in the collinear factorization framework, identifying kinematical regions 
where they are moderate or small. There one may use formulations incorporating power 
corrections from transverse quark momentum with greater confidence. This is the aim of 
the present contribution. 

In the following we show results for the convolution of the unpolarized quark and gluon 
GPDs and with the corresponding hard-scattering kernels and the asymptotic form of 
the p distribution amplitude. We model the GPDs using a standard ansatz based on double 



o 



X 



DIS 2007 



1 




Figure 1: Imaginary part of the convolution integral for the sum of gluon and quark singlet 
distributions as a function of the renormalization and factorization scale fi. 



distributions [TT], with the CTEQ6M distributions as input. Unless indicated explicitly, we 
take t = and set the factorization and renormalization scales equal, ^ ~ fip = fif(. 



2 Numerical results 
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In a wide kinematical range at small xb, 
we find huge NLO corrections which have 
opposite sign to the Born term and almost 
cancel it. This is shown for = 2 x 10^^ 
in Fig. [U where there are indications for 
an onset of pcrturbative stability at Q = 
7 GeV, but not yet at Q = 4 GeV. Taking 

= 2 X 10"'' one finds no stability even 
at Q = 7 GeV, whereas for xs = 2 x 10^^ 
the corrections arc of tolerable size already 
at Q = 4GeV. 

Figure [2] shows that in kinematics rele- 
vant for HERA measurements, NLO correc- 
tions have a huge effect on the cross section 
and moreover lead to a fiat energy behav- 
ior in conflict with experiment. Due to the 
strong cancellations between LO and NLO 
terms, the dependence on factorization and 
renormalization scale is not reduced when 
going to NLO. 

As already observed in [B] the large size 
of NLO corrections at small xb can be 

traced back to BFKL-type logarithms appearing first at NLO for vector meson produc- 
tion. Such logarithms are present in many processes (including DIS) but have a rather large 
numerical prefactor in the present case. It is to be hoped that all-order resummation of 



Figure 2: Cross section for 'y*p — > pp with 
a longitudinal photon. Bands correspond to 
the range Q/2 < fi < 2Q and solid lines 
to fi = Q. We also show the power- law be- 
havior a cx iy°-®^ (with arbitrary normaliza- 
tion) obtained from a fit to data in the range 
0.001 <a;s< 0.005 [12]. 
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Figure 3: Real and imaginary part of the convolution integral for the sum of gluon and 
quark singlet distributions for fip = ^ir = Q = 2 GeV. 
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Figure 4: Renormalization scale dependence of the real part of the convolution integrals for 
the sum of gluon and quark singlet distributions (left) and for for the difference of u and d 
quark distributions (right). 



these logarithms in the hard-scattering kernel will give perturbative stability at small xb- 

In the Xb range relevant for experiments at COMPASS, HERMES, and JLAB, we gen- 
erally find corrections which are sizable but not huge. An exception is the real part of the 
gluon and quark singlet amplitudes, where corrections become large for decreasing xb, as is 
seen in the left panel of Fig. \3\ 

In the quark nonsingiet sector there are large terms in the NLO kernel due to gluon 
self-energy corrections. The BLM procedure for setting the renormalization scale aims at 
resumming these to all orders in a^. Applied to the process at hand, one finds however that 
this requires fj.R to be typically an order of magnitude smaller than Q |131ll4j. This is outside 
the validity of the perturbative calculation for most practically relevant Q. Numerically 
we find that for fiR < 2 GeV the NLO corrections become unstable for several convolution 
integrals, as shown for examples in Fig. |4l 

We have therefore omitted this region when estimating the scale setting error in Fig. [5l 
where we show the cross section in typical fixed-target kinematics. We see that NLO correc- 
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Figure 5: Cross section for — > pp with a longitudinal photon. Bands correspond to the 
range 2 GeV < ^ < 4 GeV in the left and to 2 GeV < /i < 6 GeV in the right plot, and solid 
lines to /.( = Q in both cases. 

tions are quite large for = 4GeV^, whereas for = 9GeV^ and xb > 0.1 they become 
moderate. 
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